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Controlled non-local energy and coherence transfer enables light harvesting in photosynthesis 
and non-local logical operations in quantum computing. The most relevant mechanism of coherent 
coupling of distant qubits is coupling via the electromagnetic field. Here, we demonstrate the 
controlled coherent coupling of spatially separated excitonic qubits via the photon mode of a solid 
state microresonator. This is revealed by two-dimensional spectroscopy of the sample's coherent 
response, a sensitive and selective probe of the coherent coupling. The experimental results are 
quantitatively described by a rigorous theory of the cavity mediated coupling within a cluster of 
quantum dots excitons. Having demonstrated this mechanism, it can be used in extended coupling 
channels - sculptured, for instance, in photonic crystal cavities - to enable a long-range, non-local 
wiring up of individual emitters in solids. 



Sunlight absorption via antenna proteins and the sub- 
sequent resonant energy transfer over a few nanometers 
towards the reaction centre is at the heart of the photo- 
synthesis process 1 ' 2 . In this case the electronic or Forster 
coupling strength and the local vibrational dissipation is 
tuned to achieve a fast directional energy transfer. In- 
stead, in quantum information science and cavity quan- 
tum electrodynamics (cQED), the quantum bus tech- 
nology aims to provide dissipation-less coupling between 
distant quantum systems. Within this field, long-range 
coherent coupling between individual, distant supercon- 
ducting qubits has recently been demonstrated 3 , so that 
the construction of quantum logic gates and networks is 
within the reach of present technology 4 6 . In this con- 
text, qubits embedded in a solid state matrix are attrac- 
tive, as they can fully benefit from the lithographic and 
materials processing techniques developed in the semi- 
conductor industry to fabricate on-demand photonic lay- 
outs, so as to enable long-range coherent coupling within 
a discrete set of qubits via a photonic intra-cavity bus. 
Moreover, their "niche" with respect to superconducting 
qubits is marked out by a simultaneous functionality at 
relatively high temperatures and at optical frequencies, 
and they provide an interface to flying qubits such as 
photons in optical fibres. 

In recent years, significant progress has been made 
in the realization of high quality optical microres- 
onators which enabled pioneering demonstrations of the 
strong 7 10 and quantum strong coupling regime 11 ' 12 in 
solid state. Micropillar cavities (cf. Fig. la) are a model 
system for the study of strong coupling in this field. They 
consist of self-assembled InGaAs quantum dots (QD) - 
providing individual exciton states (X) of high oscilla- 
tor strength, located in the anti-node of the fundamental 
cavity mode (C). In these structures a quantum of opti- 
cal excitation coherently oscillates between the fermionic 
exciton and bosonic cavity photon state. The resulting 



eigenstates of mixed exciton and photon character form a 
Jaynes-Cummings (JC) ladder 13 with increasing number 
of photons in the cavity mode, showing a Rabi splitting of 
the rungs proportional to the root of the photon number. 

In this Letter we report on coherent measurements and 
modeling of cavity-mediated coherent coupling between 
three quantum dot excitons, moving from the JC ladder 
to the Tavis-Cummings (TC) ladder 14 . This constitutes 
a crucial step towards a quantum bus based on semicon- 
ductor photonic structures. 

Coherent photonic coupling of distant qubits is real- 
ized here by their dipole interaction with a common op- 
tical mode of a high-quality microresonator. For this 
purpose we have chosen a micropillar similar to that em- 
ployed in Ref. 12 as described in the Methods section. By 
micro-photoluminescence (/iPL) measurements as shown 
in Fig. lb, we have identified a triplet of Xs - labeled XI, 
X2 and X3, which at a temperature of 8K are slightly 
blue-shifted from the cavity. Increasing the sample tem- 
perature, the Xs (short-dashed lines) are tuned through 
C (long-dashed line) due to a reduction of the semicon- 
ductor band-gap. The data reveal three X-C avoided 
crossings at around 13 K, 21 K and 25 K, showing that 
each of them is in the strong coupling regime. 

A triple exciton - cavity system has a level scheme 
as shown in Fig. 2, more complex than the previously 
studied single exciton-cavity system 12 . It hosts four 
polaritonic transitions from the vacuum state. The 
polariton frequency tuning (solid lines in Fig. lb) as 
well as the variation of polariton linewidths (solid lines 
in Fig. lc) can be described by a coupled oscillator 
model with the Xl-C, X2-C, X3-C coupling parame- 
ters #2? #3) = (43,40,31.5) /ieV, homogeneous broad- 
enings (7x1, 7x2, 7x3, 7c) = (18, 11.5, 16, 36.5) /ieV, and 
frequency distances cjx2 — <^xi = 131 /ieV, ujx3 — <^xi = 
248 fieV. The parameters were obtained from a global 
fit of the coupled oscillator model to the detuning- 
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FIG. 1. Characterization of the investigated quantum dot - micropillar system, a) Sketch of the micropillar structure 
including the light coupling from the top facet, b) Temperature dependent photoluminescence spectral intensity under non- 
resonant excitation on a linear grey scale black (0) to white. The bare resonance energies of excitons and the cavity mode (white 
dotted and dashed lines, respectively), and the coupled polariton energies (solid lines) obtained from a Lorentzian lineshape fit 
and modelling (see supplementary material) are overlayed to the data. The corresponding average detuning 5 (see Eqn. 1) is 
shown on the upper axis scale, c) Coupled resonance linewidths (measured: symbols, modeling: lines). Colours and linestyles 
as in b). 



dependent transition energies and broadenings deter- 
mined by Lorentzian lineshape fitting of the /iPL-spectra 
in Fig. lb as described in the supplementary material 
(SM). To describe the detuning in this TC system with 
non-identical two-level systems, we introduce the average 
cavity detuning 

3 / 3 

5 = UJ C -Y1 9n^Xn j^Qn (1) 

71—1 I 71— 1 

having the values of 5 = (—124, 16, 124) fieV at the three 
anti-crossing points and shown in the top axis of Fig. lb. 
Multi-polaritonic features as reported in Fig. lb,c were 
previously observed in /iPL experiments 15,16 and indi- 
cate via an anti-crossing behaviour between individual 
excitons and the cavity mode that Xs could be coher- 
ently wired up by the cavity field. However, a direct 
measurement of coherent coupling is not afforded by /iPL 
measuring the incoherent emission, neither does it pro- 
vide means to perform coherent optical manipulation in 
prospective quantum bus nanophotonic structures. An 
explicit demonstration and selective manipulation of the 
cavity- mediated coherent coupling of excitons requires 



extracting and controlling the coherent response of the 
photonic resonator. 

To demonstrate and study the cavity mediated coher- 
ent coupling of Xs we employ here heterodyne spectral 
interferometry 17 (HSI) to detect the cavity-mediated co- 
herent coupling between the excitons. This technique 
is used to measure the four-wave mixing (FWM) of the 
strongly-coupled Xs-C system, which arises from the op- 
tical nonlinearity of the quantum dot excitons. In short, 
two pulses Ei and E2 of 1 ps duration and variable delay 
r are exciting the fundamental cavity mode, see Fig. S6 
in the SM. The resulting FWM polarization emitted from 
the micropillar P(t,r) ex E\E\ and its Fourier transform 
(FT) versus t, P(cj,t), are measured using spectral in- 
terferometry. 

The coherent dynamics giving rise to FWM involves 
the four Xs-C mixed states of the first rung of the TC 
ladder, which were identified in /iPL, and additionally 
the 7 polariton states of the second rung which are mixed 
states of the family of uncoupled exciton-photon states of 
the second rung: one two-photon state without excitons, 
three one-photon states with one of the three excitons 
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FIG. 2. Level scheme of the Tavis-Cummings ladder of the 
three exciton - one cavity system, and transitions relevant for 
the coherent FWM response, for 5 = — 29/xeV. a) Coherence 
created by the pulse arriving first {E\ for r > 0, E 2 for r < 0). 
b) Transitions emitting FWM after the arrival of the second 
pulse. 



filled, and three zero-photon states with two of the three 
exciton states filled. A sketch of the resulting ladder of 
levels at T=19K (6 = -29/ieV) is shown in Fig. 2, to- 
gether with the transitions relevant for the delay-time 
and real-time coherent dynamics probed in FWM. For 
positive delay r > 0, pulse E\ arrives first and creates 
a one-photon coherence given by a wavepacket of the 4 
states of the first rung, which is coherently evolving until 
the arrival of E 2 . Conversely, for negative delay r < 0, 
pulse E 2 arrives first and creates a two-photon coher- 
ence given by a wavepacket in the second rung, which 
is coherently evolving until the arrival of E\. In both 
cases, at the arrival of the second pulse FWM is created 
as a superposition of all optical transitions between the 
ground state, first and second rung, consisting of 4 tran- 
sitions between the ground state and the first rung and 
28 transitions between the first and second rung of the 
TC ladder. We calculate the FWM polarization analyt- 
ically taking into account the states up to the second 
rung, by solving the master equation for the density ma- 
trix using a standard Xs-C coupling Hamiltonian and a 
Lindblad dissipation operator (see SM). This approach 
is exact for the third-order FWM signal of an initially 
unexcited system. 

The resulting measurements and corresponding pre- 
dictions of FWM of the system are given in Fig. 3 as 
function of time delay r for two different detuning pa- 
rameters. For S = —29/ieV (T = 19 K), the spectrally 
resolved \P(cj, r)| 2 and time-resolved |P(t, r)| 2 are shown 
in Figs. 3a,c, respectively. A dynamics significantly richer 
than in a single exciton case 12 is observed, as expected 
from the larger number of levels in the first and second 
rungs, providing 32 instead of 6 transitions contribut- 



ing to the FWM (see Fig. 2). The time-integrated FWM 
power \P\f nt (r) = j\P(t,r)\ 2 dt and the power \P(t m , r)| 2 
at a given time t m = 21 ps corresponding to the build-up 
lag of the FWM in such strongly coupled exciton- cavity 
systems 12 are presented in Fig. 3e. On a qualitative level, 
we notice the FWM beat as a function of r with a pe- 
riod of about 17 ps, corresponding to a spectral splitting 
of 243 /ieV. This is much larger than the Rabi splitting 
of any individual X, and is close to the total splitting of 
+ 92 + 93) = 229 /ieV, indicating that all four po- 
laritons contribute towards the coherent dynamics. In 
Figs. 3b, d we present the predicted FWM corresponding 
to 3a, c, using the exciton and cavity parameters retrieved 
from the jaPL data (see Fig. 1 and the SM). The predic- 
tion, which takes into account the coherent evolution in 
the TC ladder shown in Fig. 2, reproduces the rich fea- 
tures of the measurements quantitatively. 

Modifying the detuning, we can adjust the system to 
exhibit only one exciton in resonance with the cavity, 
while the other excitons are significantly detuned, so 
that the dynamics resembles that of a simpler single- 
exciton cavity system. This is achieved at S = —133 fieV 
(T=13.5K), for which XI is in resonance with C within 
5 /ieV, whilst X2 and X3 are detuned by 135 /ieV and 
253 /ieV, respectively. The measured and predicted 
\P(uj,t)\ 2 and |P(t,r)| 2 are presented in Fig. 3f,h and 
Fig. 3g,i, respectively, while the time-integrated FWM is 
displayed in Fig. 3j. We observe a beat versus delay r 
with a period of about 50 ps, corresponding to a polari- 
tonic splitting of 83 /ieV, somewhat larger than the calcu- 
lated splitting of 60 /xeV, which is slightly below 2gi due 
to the finite damping. The faster than expected beat pe- 
riod is due to remaining influence of the two additional 
excitons, as shown by the agreement of the predicted 
dynamics including all excitons (see solid line) with the 
measurements. 

The agreement between the measured FWM dynamics 
and independently predicted FWM in the framework of 
the cavity-mediated coupling model for different detun- 
ing is revealing the coherent wiring up of the three Xs via 
the cavity mode in the studied microresonator. A defi- 
nite display of the coherent coupling is afforded by the 
two-dimensional (2D) frequency domain representation 
of the FWM, in which coherent coupling is observed as 
off-diagonal signals 1 ' 2 ' 18 . We retrieve the 2D FWM 17 by 
Fourier-transforming P(lu,t) from the delay time r into 
the conjugated frequency uj t yielding P(uj,uj t ). In this 
transformation we use only positive delays r > 0, such 
that uo T represents the frequency of the first-order po- 
larization created by E\. The resulting two-dimensional 
FWM diagrams P(uj,uj t ) are presented in Fig. 4a, d. In 
this representation, coupled resonances manifest them- 
selves by corresponding off-diagonal components, show- 
ing that resonances excited by the first pulse in first-order 
(uj t axis) coherently couple to a different resonances emit- 
ting in third-order (uj axis) . To enable the Fourier trans- 
formation we set the relative phase of data at different 
delays r using a phase correction 19 at a given u (see yel- 
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FIG. 3. Delay time dependence of coherent response for T = 19 K (top) and T = 13.5 K (bottom). Spectrally resolved FWM 
power |P(cj,t)| 2 , measured (a,f) and predicted (b,g), on a logarithmic colour scale over 4 orders of magnitude. Time-resolved 
FWM power |P(£,t)| 2 , measured (c,h) and predicted (d,i), over 3 orders. (e,j) Time-integrated FWM power |P| 2 nt( r ) 5 measured 
(black circles) and predicted (red line), and measured |P(21ps,r)| 2 (blue triangles). The noise of |P| 2 n t(^) is given as open 
circles. 



low arrows) as discussed in the SM. The effect of the 
phase-correction is illustrated in the difference between 
the predictions with and without phase correction given 
in Fig. 4b, c, respectively. The phase correction leads to a 
weighting in the first-order frequency, but does not qual- 
itatively change the off-diagonal structure. 

Due to the contributions of the 28 second-rung tran- 
sitions, the features in 2D FWM are broadened. Post- 
selecting the FWM signal emitted after a survival time 
t s after its creation by the second pulse at t = 0, resulting 
in the 2D FWM P s (uj, u; r , t s ), we can suppress the fast de- 
caying components 12 . The measured and predicted post- 
selected P s (u;,cj r ,42.5ps) (see Fig. 4d,e,f) are dominated 
by the first-rung transitions and display a clear separa- 
tion between the multiple off-diagonals at the polariton 
frequencies Ai^ (see SM), showing cavity- mediated mu- 
tual coherent coupling between Xs. The corresponding 
data for T = 13.5 K (6 = -133 fieV) given in the SM, 
Fig. S6, shows the coherent coupling between two polari- 
tonic modes dominated by the cavity and XI. A de- 
tailed analysis of the strength of the off-diagonal peaks 
and their relation to the coherent coupling strength will 
be presented in a forthcoming work. 

The observed coherent coupling between three qubits 
via the cavity mode demonstrates that an optical mi- 



crocavity can act as a coupling bus for excitonic qubits. 
The coherent interaction can be controlled by non-local 
tuning of the cavity mode or the exciton energy, or by 
switching the excitonic qubit between bright and dark 
exciton states 20 . The coupling via the cavity mode 
has the prospect to be spatially extended in coupled 
cavity structures, for example in photonic crystals 21 . 
The small size of the quantum dots with respect to 
the cavity enables coupling of many more than three 
qubits, enhancing the prospects of quantum information 
processing on a chip with a photonic interface. 

Methods 

Samples: The fabrication process of high-Q quantum 
dot micropillar structures starts with the growth of 
a planar microcavity on an undoped GaAs substrate 
by molecular beam epitaxy. A GaAs A-cavity with a 
single layer of self-assembled Ino.4Gao.6As QDs with 
a ground state exciton emission wavelength around 
930 nm is sandwiched between two distributed Bragg 
reflectors (DBR) providing a vertical photon confine- 
ment. The upper (lower) DBR consists of 26 (30) 
mirror pairs, which are made of A/4 layers of GaAs 
and AlAs. For lateral photon confinement micropillar s 
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FIG. 4. Two-dimensional FWM at T=19K. Power 

|P(o;, cj t )| 2 , measured (a) and predicted with (b) and with- 
out (c) phase correction. Logarithmic colour scale from 0.09 
(black) to 1 arbitrary units. Yellow arrows indicate the uj of 
the phase correction. Magenta ticks indicate the polariton 
frequencies of the first rung Ai,fc (see SM). (d,e,f) as (a,b,c), 
but showing the post-selected |P s (cj, c<; r , 42.5ps)| 2 on a linear 
colour scale from 0.01 to 0.24 arbitrary units. 

with a nominally circular cross-section and diameters 
between 1.6 /im and 2 /im were structured by means of 
high resolution electron beam lithography and electron 
cyclotron etching. The fundamental cavity mode in the 
investigated structure has a Q-factor of 30000, and a 
corresponding photon lifetime of about 10 ps within its 
mode volume of approximately 0.3 /im 3 . 

FWM experiment: We employ heterodyne spec- 
tral interferometry 17 , using a pair of pulse trains, 
Ei and E<i at a repetition rate of 76 MHz and 1 ps 
pulse-duration which are frequency up-shifted by 
(£V2tt, (V2tt) = (79, 80.7) MHz by acousto-optic 
devices. The spatial mode of E\ and E2 is matched to 
the cavity mode on the top facet of the micropillar by a 



microscope objective mounted inside a helium bath cryo- 
stat. The pulses are resonant to the cavity mode and are 
injecting into the intra-cavity field in average 0.25 (0.75) 
photons per pulse of E\ (E2), respectively, see Fig. S7 
in SM. The emitted light is interfered with a frequency 
unshifted reference field and the FWM signal is selected 
at the heterodyne beat note 2Q 2 — ^1 and spectrally 
resolved. The complex FWM response P(t,r) oc E\E\ 
is retrieved by spectral interferometry 17,22 . 

Theory: We use a standard approach 15,16 ' 23 , with 
the Hamiltonian of a single photonic mode coupled to 
three two-level systems with dissipation of all compo- 
nents taken into account by a Lindblad super-operator. 
To calculate the FWM polarization, we follow our earlier 
rigorous approach 12,24 to the optical response of a system 
on a sequence of optical pulses. For the present system 
excited by two ultrashort pulses it has the analytic form 

PfwmM = ^2c jk e^ t e iXkT (2) 

jk 

in which = AJ fc (A& = \2k) are the complex polariton 
frequencies of the first (second) rung, for positive (nega- 
tive) delay times r, respectively, while ojj are the complex 
frequencies of all possible transitions between the ground 
state and the first rung and between the first and second 
rung, all contributing to the real-time dynamics. 

The expansion coefficients Cjk are calculated exactly 
reducing the full master equation for the density matrix 
to a finite matrix problem. Further details are available 
in the SM. 
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I. FOUR- WAVE MIXING RESPONSE OF N QUANTUM DOTS COUPLED TO THE 

CAVITY 

The Tavis-Cummings (TC) model of N quantum dots (QDs) coupled to a single photonic mode in a 
microcavity is described by the Hamiltonian 



N 



n=l 



uxn\n)(n\ +^ n (a t |0)(n| +a|n)(0|) 



(1) 



in which ujq is the frequency of the cavity mode, uoy^ is exciton transition energy of the n-th QD, and 
g n is the coupling strength between the n-th QD excition and the cavity mode. Furthermore a) (a) is 
the cavity photon creation (destruction) operator, |0) is the ground state (GS) of the exciton system, 
and |n) denotes the state with one exciton in the n-th QD and all other QDs empty. We use ft = 1 for 
simplicity of notations. We assume that all QD excitons are close to resonance with the cavity mode, 
and neglect any excited or multi-excitonic states within the same QD, as they are generally off resonance. 
However, since we consider QDs which are electronically uncoupled, multi-excitonic states in which all 
excitons are in different QDs stay close to resonance with the cavity mode, and significantly contribute 
to the four- wave mixing (FWM) dynamics. 

The system of the N QD excitons coupled to the cavity is excited by an external electric field consisting 
of two ultrashort pulses separated by the delay time r and having complex pulse areas \iE\ and pE 2 , 
where /i is the effective dipole moment of the cavity mode. The measured FWM polarization is a third- 
order signal proportional to E\E\. In the calculation, the effect of the pulses is taken into account as 
follows. The first pulse E\ contributes in first-order, changing the density matrix according to 

pW=-ipE*[a,p^}, (2) 

where p^ (p^) is the density matrix before (after) the pulse. The second pulse E 2 contributing in 
second order results in the following change of the density matrix: 

pW^^hM^rt. (3) 

Between and after the pulses the time evolution of the density matrix p(t) is described by the master 
equation 

(4) 

where L is the Lindblad super-operator of the exciton-cavity system, 

N 



Lp = [H, p] - Z7c (a^ap + pa) a - 2apo)) + ^ 7xn 



n=l 



\n)(n\p + p\n)(n\-2\0)(n\p\n)(0\ 



(5) 



with 7c and 7xn being, respectively, the cavity damping and the exciton dephasing, which are assumed 
to be lifetime limited. Equation (4) has the following formal solution 

p(t)=e- iL ^p{t>), (6) 
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and thus the FWM polarization takes the form 



P(t,r) = Tv{pa} = Li^V^Ef • Tv{R(t, r)a} 



where the third-order component of the density matrix is proportional to 



R(t,r) 



iLt 



a\[a\e- itT [a,p(-oo)]] 



for positive delay between the pulses (r > 0), and to 



R(t,r) 



-iLt 



a,e iLr [a^[a^p(-oo)]\ 



(7) 



(8) 



(9) 



for negative delay (r < 0). Here p{— oo) is the density matrix before the excitation. 

We solve equations (2)-(4) exactly by expanding all operators into a set of uncoupled exciton-photon 
states \i) = |nxi, n x2 , . . . , nx/v; nc) (i = 0, 1 , £ . . . ), where nxn = 0, 1 and nc = 0, 1, 2 . . . are, respec- 
tively, the exciton and photon occupation numbers. We assume that the system is initially in its ground 
state, with the density matrix p(— oo) = \ 0)(0\. This is justified by the excitations of the system having 
an energy about 3 orders of magnitude larger than the thermal energy and coherence times shorter than 
the period of the repetitive excitation (13 ns) used in this work. The resulting FWM dynamics includes 
only transitions between the GS, the first rung, and the second rung of the TC ladder. With TV quantum 
dots coupled to the cavity, there are N± = 1 + N states in the first rung and 7V 2 = 1 + N(N + l)/2 states 
in the second rung, so that the expansion of density matrix has the form 



Pit) 



N!+N 2 



(10) 



i,j=Q 



with 1 + Ni + N2 basis states. A similar expansion is used for the photon creation (annihilation) operator 
contributing to Eqs. (2) and (3). To solve the master equation (4) after the pulsed excitation, a bigger, 
M x M matrix of the super-operator L is introduced, with M = Ni(l + N2) being the total number of 
all possible transitions, between the GS and the first rung (iVi), and between the first and second rungs 
(N1N2). Only such transitions contribute to the FWM polarization. The solution to Eq. (4) then has the 
matrix form 



pit) = e~ tLt p(0) = Ue- %ut Vp(0) , 



(11) 



in which the matrix L is diagonalised as 



l = unv , UV = 1 . 



(12) 



with U and V being matrices of right and left eigenvectors and Cl a diagonal matrix of complex eigenfre- 
quencies ujj of all possible transitions involved. The density matrix then becomes a finite superposition 
of exponentials, 



M 



p(t) = E 



r 7 e 



with the frequencies Cjj and matrices rj calculated exactly. In particular, 



Xi^k and ujj = A 2 ,fc - AJ >r: 



(13) 



(14) 



for the GS -first rung and first rung-second rung transitions, respectively, where Ai^ (^2,k) are the 
complex energy levels of the first (second) rung, calculated by diagonalizing of the complex symmetric 
matrix of a non-Hermitian Hamiltonian (see below). The time evolution between pulses is given by a 
smaller matrix L as only transitions from the GS to the first (second) rung participate in the delay 
dynamics, for positive (negative) delay times. The FWM polarization then takes the following explicit 
form: 



M 



k=l 
N 2 



(15) 



k=i 
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With the analytic form Eq. (15) at hand, it is straightforward to calculate the Fourier transform (FT) of 
the FWM polarization. In particular, for positive delays, its two-dimensional FT has the form 



M iVi 



a-ik 



(16) 



where a Lorentzian spectrometer resolution 75 (half width at half maximum) is included. 

Although our analytical approach is general and can be used for any arbitrary N producing a reasonable 
size of the basis (the size of the matrix L scales as A 7 " 3 ), in the present calculation we concentrate on a 
system of TV = 3 QDs coupled to the cavity. Then the basis reduces to the following twelve states which 
include Ni = 4 states of first rung and N 2 = 7 states of the second rung: 

Ground state |0) = |0,0,0;0); 

First rung |i ) = |0, 0, 0; 1} , \2) = |1, 0, 0; 0) , \S) = |0, 1, 0; 0) , \4) = |0, 0, 1; 0); 



Second rung 1 5) = |0, 0,0; 2), \6) = |1, 0, 0; 1) , 1 7) = |0, 1, 0; 1) , \8) = |0, 0, 1; 1) , 

|S) = |1,1,0;0), |i0) = |l,0,l;0), \11) = |0, 1, 1; 0) . (17) 

Both the density matrix and the photon creation operator are 12 x 12 matrices, while the super-operator L 
after the excitation by both pulses is represented by a 32 x 32 matrix which is due to M = 32 transitions: 
4 transitions between the GS and the first rung and 4 x 7 = 28 transitions between the first and second 
rungs. 

The energy levels of the TC ladder are calculated by diagonalizing the effective non-Hermitian 
Hamiltonian: 



H = 



/ffi ...\ 

H 2 .., 

H 3 ... 

V ! ! ! 'J 



(18) 



The first two blocks on the main diagonal of H refer to the first and second rungs, respectively: 
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where cD c = ^c ~ he and cDxn ^Xn - ^7xn- 



II. PHOTOLUMINESCENCE CHARACTERIZATION 



The investigated sample consists of fields of 6 by 6 micropillars of nominally equal diameter, with 50 /im 
separation in x and y direction corresponding to the [110] and [110] crystallographic directions of the 
GaAs substrate, respectively. Each micropillar was accompanied by a reference pillar of 10 jam diameter 
at 10 /im distance which was used to reflect the reference beam in the heterodyne spectral interferometry 
(HSI) experiment. A scanning electron microscopy image of a row is shown in Fig. Sla. 

The micropillars have diameters of 1.5, 1.6, 1.7, 1.8 and 1.9 /im, and for each diameter 10 equal fields 
provided 360 nominally equal micropillars. The sample was mounted on an x — y stage in a helium 
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SEM optical 




FIG. SI. Images of the investigated sample, a) Scanning electron microscopy image of a row of micropillar - 
reference pillar pairs, b) White light reflection image taken in the cryostat at 10 K using the same optical setup 
as in the confocal PL and HSI experiments. Scale bars are given. 



bath cryostat, together with a microscope objective of 0.85 numerical aperture which was adjustable in 
all three dimensions by a home-build stage using piezoelectric bender actuators. A reflection image of 
a micropillar taken with this arrangement is shown in Fig. Sib. To identify micropillars in the strong 
coupling regime, we used confocal photoluminescence (//PL) excited by a laser of 532 nm wavelength, 
focussed to a spot of about 0.5 fim FWHM on the top face of the micropillar. Typically 10% of the 
micropillars with 1.5 — 1.7 /im diameter showed strong coupling. The micropillar presented in this work 
was selected from about 200 micropillars for the strongest coupling of multiple quantum dots. The 
photoluminescence excitation power at the micropillar was about 100 nW. At lower intensities the PL 
spectrum remained essentially unchanged, while at higher intensities line broadening and reduction of 
the Rabi-splitting was observed. 

To characterize the quantum dot - cavity coupling, /iPL spectra were measured as function of the 
sample temperature from 7K to 35 K. The PL spectrum at each temperature was fitted by a sum of 
Lorentzian lines, as exemplified in Fig. S2. The resulting line positions as function of temperature were 
fitted with the first rung transition energies Ai^ given by Eq.(19). In the fit, we used an explicit functional 
dependence for the temperature-tuning of the energies of quantum dots and cavity, from a model of the 
band-gap shift of semiconductors 1 . 

uxn(T) = uxniO) + F(T), w c (T) = w c (0) + ^(T) (21) 

with 

^>=-tM^H (22) 

The parameters describing the temperature-dependence were determined from a fit to the emission 
of quantum dots detuned from the cavity by more than 1 meV, yielding a — (60.9 ±0.6) /ieV/K and 
6 = (58.9 ± 1.0) K. The quoted errors represent the statistical variation of the results between different 
fitted quantum dots. 

From the fit of Ai^ to the measured line positions and widths, we deduce the low temperature energies 
^Xn(O), cjc(0) 5 the linewidths 7xn?7c and the coupling strengths g n , as shown in table I. Additionally, 
the scaling factor r] = 0.227 ± 0.003, for the cavity shift was determined, showing that the cavity shift is 
about 4 times less than the quantum dot exciton shift. 

For simplicity, we assume here that the linewidths and the coupling strengths are independent of 
temperature, and that the linewidths are Lorentzian, i.e. homogeneously broadened. A spectrometer 
resolution of 75 = 4 fieV (HWHM) was subtracted from the fitted linewidths. 

The intrinsic homogeneous linewidth of the quantum dots 2 is expected to be limited by radiative 
decay and phonon scattering. The radiative decay time of the quantum dots in bulk GaAs is about 
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1.3342 1.3344 1.3346 1.3348 1.3350 

photon energy (eV) 

FIG. S2. Photoluminescence spectra from the investigated micropillar at a temperature of 13.5 K (top) and 19 K 
(bottom), and fits by multiple Lorentzian lines (solid lines). 



parameter 


value 


unit 


(0) 


1.3346106 


eV 




1.3347412 


eV 




1.3348584 


eV 


wc(0) 


1.3345732 


eV 


9i 


43 


fieV 


92 


40 


fieV 


93 


31.5 


fieV 


7xi 


18 


fieV 


7X2 


11.5 


fieV 


7X3 


16 


fieV 


7c 


36.5 





TABLE I. Parameters of the three-exciton cavity system deduced from a fit to the temperature-dependent pho- 
toluminescence spectra. 



400 ps, corresponding to 1 fieV HWHM. In micropillars, this decay is not significantly modified due to 
the presence of leaky modes which have a similar local density of states as the bulk modes 3 . The phonon- 
scattering increases the zero-phonon linewidth for quantum dots with 140 meV confinement energy 2 by 
about 0.5 /xeV at 20 K and 2/ieV at 30 K. 

The measured linewidths of different uncoupled quantum dots are between 10 and 30 /xeV, and vary 
from dot to dot, both in terms of low-temperature linewidth and dependence on temperature, which we 
attribute to spectral diffusion. Such a broadening mechanism is presently not taken into account in the 
FWM modeling, and might explain some of the remaining differences between the experimental data and 
the modeling, specifically for negative delay due to the expected photon echo formation. We plan to 
extend the modeling to include also inhomogeneous broadening and pure dephasing of the quantum dot 
excitons 4,5 in future work. 

The coupling strength is expected to decrease with temperature proportional to the square root of 
the zero-phonon line weight 2 , resulting in a reduction of about 5% from 8K to 30 K. We have neglected 
this dependence. The resulting fit (see Fig. 1) reproduces the energy positions and linewidths to within 
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the measurement error of about 5 jaeV. The above approximations are therefore resonable for the data 
at hand. The linewidths of the strongly coupled modes are actually dominated by the cavity linewidth, 
which is given by the photon lifetime and thus Lorentzian. This reduces the importance of the broadening 
mechanism in the quantum dot excitons. 



III. POLARITON STATES 
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FIG. S3. Energies and linewidths of the polariton states of the first and second rungs and their probability distri- 
bution (grey scale matrices) in the basis of the uncoupled states Eq. (17) of the exciton-cavity system, calculated 
for (a) model symmetric structure and (b) realistic asymmetric structure with S = — 29 /xeV (corresponding to 
T — 19 K). The parameters of the system used for the symmetric structure are g n = 40 /xeV, 7xn = 0, 7c = 60 /xeV, 
cjx2 = cjc, and cjx2 — cjxi = cjx3 — cjx2 = 130 /xeV (corresponding to zero detuning). The parameters of the 
realistic structure are given in the main text and table I. Here, the similar contributions of the cavity mode to all 
states results in similar linewidths of the levels of a given rung. 

The polariton states in the micropillar calculated for a symmetric model structure and realistic param- 
eters of the investigated sample are shown in Fig. S3. The detuning dependence of the cavity and exciton 
levels are shown in Fig. S4. 

IV. SUPPLEMENTARY FWM RESULTS 

Here we show supporting experimental and theoretical results concerning the four- wave mixing. Adding 
to the data shown in Figs. 3,4 of the main manuscript, we show here the calculated FWM dynamics for 
various detunings in Fig. S5, both the real-time resolved |P(t,r)| 2 , the frequency resolved |P(cj,t)| 2 , and 
the 2D FWM \P(u, u T )\ 2 . 

The phase correction applied to the experimental FWM polarization P(uj,t) setting the phase evolution 
of the FWM versus r to its expected evolution for an uncoupled system, is given by 



P cor (^,T) = P(cj,r)exp (i (u COT r - arg P(u cor , r)^ . 



(23) 
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FIG. S4. Energies of the uncoupled QD-exciton and cavity modes as function of the detuning S = ujq — 
(Xn=i dn^Xn) / ($^n=i 9 n ) ■ Dotted vertical lines correspond to the values of the detuning used in Fig. S5. 

Coherent coupling analysis for 5 = —133 jaeV (T = 13.5 K) is shown in Fig. S6. (a) 2D FWM diagram 
|P cor (u;, uo T )\ 2 obtained by phase correcting P(lj,t) at u cor = 1334.52 meV (as marked by the arrow) and 
Fourier-transforming the corrected data along r. Dash-dotted line represents the diagonal cj = uo T . The 
2D FWM response contains both first- and second-rung signals masking the coherent coupling features 
due to the broad spectral width of second- to first-rung transitions. Theoretical spectra corresponding 
to (a) are calculated with (b) and without (c) phase correction. In (a-c) linear colour scale over 2 orders 
of magnitude is used, (d) 2D FWM retrieved by using the signal emitted after a time-lag of 42.5 ps, and 
the phase correction routine is applied as in (a). In this case the 28 second- to first-rung transitions are 
virtually suppressed, so that the FWM signal is dominated by the 4 first-rung to ground state transitions. 
An off-diagonal resonance caused by the cavity mediated coherent coupling of the QDs is detected and 
well reproduced by the corresponding theoretical simulation, with (e) and without (f) phase correction. 

To investigate the nonlinear regime of our measurements, we have performed FWM measurements as 
function of the excitation power in the pulses E\ and E 2 at zero delay time r. The resulting FWM 
spectra and integrated power are shown in Fig. S7. A third order scaling is observed, as expected for the 
lowest order FWM signal, which is proportional to E\E\. For higher excitation, the third order scaling 
saturates, followed by a strong reduction. The total driving power used in the experiment (marked by 
the vertical line) is within the third-order regime. This implies that polariton levels involved in the FWM 
experiment are limited to only the first and the second rungs of the TC ladder. We estimate that in 
the conditions of our experiment each E\ (E 2 ) pulse delivers 80 (240) photons onto the top facet of the 
micropillar. This corresponds, in average, to 0.25 (0.75) photons per pulse that are injected into the 
cavity. With increasing power, a transition to a response dominated by the cavity resonance is observed, 
resembling a Mollow triplet. We observed an evolution towards a Mollow triplet in a single X-C system 6 . 
Modeling of the high-excitation response will be reported in a future work. 
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FIG. S5. Calculated FWM versus detuning using the parameters listed in table I. FWM intensities 
\P(t,r)\ 2 (left column), \P(oj,t)\ 2 (middle column), and \P(oj,oo t )\ 2 (ri ght column), for the detunings 5 = 
(0, — 40, — 80, — 120)/ieV from top to bottom. Logarithmic colour scale as shown over 4 orders of magnitude. 
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FIG. S6. 2D FWM at T = 13. 5K (8 = -133/xeV). (a) \P(u, cj t )| 2 , phase corrected at cj cor = 1334.52meV, 
marked by the arrow. The diagonal uo — uu T is shown as dashed line. Predicted data corresponding to (a) are 
shown with (b) and without (c) phase correction. In (a-c) a linear colour scale over 2 orders of magnitude is used, 
(d) post-selected |P S (<^, cj t , 42.5 ps)| 2 , and its prediction with (e) and without (f) phase correction. 



10 



# of photons per pulse E 1 entering the cavity 

0J 1 

— i — i i 1 1 1 1 1 1 — i — i i 1 1 1 1 1 1 — i— 



reaching the top facet of the sample 



10 



100 



1000 -C 



10 b - 



1 10* 

3 



10 d r 



CD 

I 10 2 

Q. 



^ 10' 

LL 

10° 



irq 1 1 — i i i i 


rq 1 1 — i i i i 1 1 1_ 




• 


■ 




: 


• 


r 

• 


- 

T=19K | 




I^si^i 2 ] 


• 




• 




: ni 


>i i : 


0.1 


1 10 



power (|nW) 



1 — i — r 
hd 



-> 1 ■ 1 ■ r 

b A 4.7uW 





1 r 

l-e 




A 



2.46uW- 



-> r 

1.01uW- 



0.51 uW 



0.25uW" 



g 




0.14p.W. 




1334.4 1334.6 1334.8 

photon energy (meV) 



FIG. S7. Intensity dependence of the FWM power, (a) Spectrally integrated FWM power as a function of the 
combined E\ and E 2 driving power, at r = and 8 — — 29/ieV (T=19K). (b-g) Spectrally resolved FWM power 
for different excitation power as labeled. 



